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ABSTRACT 

We present the optical/infra-red lightcurve (0/IR) of the black hole X-ray 
binary GX 339-4 collected at the SMARTS 1.3m telescope from 2002 to 2010. 
During this time the source has undergone numerous state transitions includ- 
ing hard-to-soft state transitions when we see large changes in the near-IR flux 
accompanied by modest changes in optical flux, and three rebrightening events 
in 2003, 2005 and 2007 after GX 339-4 transitioned from the soft state to the 
hard. All but one outburst show similar behavior in the X-ray hardness-intensity 
diagram. We show that the 0/IR colors follow two distinct tracks that reflect 
either the hard or soft X-ray state of the source. Thus, either of these two X-ray 
states can be inferred from 0/IR observations alone. From these correlations 
we have constructed spectral energy distributions of the soft and hard states. 
During the hard state, the near-IR data have the same spectral slope as simul- 
taneous radio data when GX 339-4 was in a bright optical state, implying that 
the near-IR is dominated by a non-thermal source, most likely originating from 
jets. Non-thermal emission dominates the near-IR bands during the hard state 
at all but the faintest optical states, and the fraction of non-thermal emission 
increases with increasing optical brightness. The spectral slope of the optical 
bands indicate that a heated thermal source is present during both the soft and 
hard X-ray states, even when GX 339-4 is at its faintest optical state. We have 
conducted a timing analysis of the light curve for the hard and soft states and 
flnd no evidence of a characteristic timescale within the range of 4-230 days. 

Subject headings: binaries: close — black hole physics — infrared: stars — radiation 
mechanisms: non-thermal — radiation mechanisms: thermal — stars: individual (GX 
339-4) 
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Introduction 



Observations of X-ray emitting binary stars provide clear evidence for the presence of 
accreting black holes in some of these systems since dynamical measurements of the mass 
of the accreting compact object are greater than the maximum mass of a neutron star 
(Remillard & McClintock 2006 and references therein). Around 20 black hole X-ray binaries 
(BHXBs) have been discovered in our Galaxy so far. Observations of these objects are seen 
to have variable emission at all wavelengths due to changes in the mass accretion rate, the 
appearance/disappearance of a jet, and variations in the jet emission strength. 

GX 339-4 is unusual among BHXBs in that it has not been observed to go into an 
optical/infrared "quiescent" state during which the flux from the companion star dominates 
over the flux from the accretion flow. Observations of the optical counterpart during its 



faintest optical and X-ray state ( IShahbaz. Fender, fc Charles 



200ll ) do not show absorption 



lines in the optical spectrum, and identification of the companion star has remained elusive. 



This has complicated e 



(e.g., 



forts to determine the mass and binary parameters of the system 



Buxton &: Vennes 



20031). However, observations of emission features associated with 



the irradiated companion star and accre tion disk were succ essful in determining the orbital 



20031 ) so it is clear that this source 



period and mass function of the source (JHynes et al. 
does contain a black hole. 

While the lack of a quiescent state has complicated studies of the binary parameters 
of GX 339-4, its continuous accretion has made it a key source for studying accretion 
phenomena. To do so requires simultaneous multiwavelength data that covers the spectral 
energy distribution (SED) as broadly as possible. Th e optical and near- infrared (0/IR) 



spectral region is particu l arly important in GX 339-4 ( ICorbel fc Fender 



2006 



Homan et al. 



2005 : 



Coriat et al. 



2002; 



Russell et al. 



20091 ) as it appears to be the location of a cross-over 



between a non-thermal and thermal dominated spectrum. This is also true for other 
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BHXBs including XTE J1550-564 flJain et al 



Chatv. Dubus. &: Raichoor 



Brocksopp et al. 



2011 



Russel 



2010D and 4U 1543-47 



et al. 



20011: ICorbel et all 120021 : iRussell et al.l 12010 : 



2011 



LXT 



Buxton fc Bailyn 



J1118+480 flHvnes et al. 



2006 



2004h . 



During the soft-state of GX 339-4, the spectral slope of the 0/IR emission is positive 
suggesting that its origin is thermal in nature. The 0/IR spectrum in BHXRBs during this 
state is expected to be dominated by a multicolor blackbody representing the accretion disk. 
Such a spectrum consists of a Rayleigh- Jeans tail (a = 2, where F^ oc z/"), a flattened profile 
(a ~ 1/3 for a viscous heated disk), followed by the Wein tail. To discern the temperature 
of the thermal source one needs to observe the turnover to the Wein tail. This turnover has 
not be en directly observed in G X 339-4 but simultaneous 0/IR and UV observations of GX 



339-4 f CadoUe Be 



of other BHXBs (jHynes et al 



et al. 



2011 



) show the flux increases from the infrared to the UV. Studies 



2002 



Hynes 



2005 



Zurita Heras et al 



201 ll ) have revealed 



that the broadband SED increases from optical to UV wavelengths and, in some cases, 
seen to peak in the UV. However, the UV flux determination is sensitive to the reddening 



estimate which can be frought with difficulties ( IFitzpatrick 



19991 ) placing large errors on the 



peak location. The optical spe ctral s 



observed to be 0.5 < « < 1.5 Mynes 



ope in many BHXBs during X-ray outbursts has been 



20051 ). However, the optical spectra were insufficient 
to determine whether the disk flux arose from irradiation from an X-ray source (central or 
vertically extended) or viscous heating in the outer disk. For this one needs simultaneous 
UV data. Not surprisingly, therefore, there is conflicting e vidence as to whether the disk in 
GX 339-4 is irradiated or not during X-ray outbursts (e.g. 



2005 



Markoff et al. 



2005 : 



Maitra et al. 



2009 



Hviies et al. 



CadoUe Bel et al. 



201lh . 



2003 



Homan et al. 



In the hard state, the near-infrared spectr al slope in GX 339 -4 differs significantly 



from the optical and often 



though not always, 



spectrum (I Corbel &: Fender 



2002 



Homan et al. 



Coriat 



2OO5I : 



et al. 



2009h 



Russell et al. 



i nes up with the ra dio 



2006 



Coriat et al 
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20091 ) indicating that the majority o: 



be non-thermal emission from a jet (JFender 



near-infrared emission during the hard sta t e may 



2001 



Fender. Belloni. fc Gallo 



2004 



20051). 



The location of the non-th ermal power- la w break is aii impor tant parameter as it is used 



to estimate the jet power (JFender 



20061). 



MotchetaL 



(119851 ) first showed that, during 



the hard state, the IR to X-ray SE P was well explai n ed by a single power-law fit with a 
slope of a = —0.58. Subsequently, iCorbel fc Fended (120021 ) constructed two SEDs with 



q uasi-simult aneous radio . O / IR and X-ray data and suggested that the power-law breaks in 



the near-IR. 



Nowaketal 



( 120051 ) found that, using 10 simultaneous radio/X-ray hard-state 
observations of GX 339-4, the hard-state SED could be fit with a doubly-broken power-law 
that breaks in the IR, and that the position of the break changes with X-ray flux. Recent 
studies are now sho wing this to be the case. Using simultaneous radio and SMARTS 0/IR 



data taken in 2010, 



Coriat et al. 



( 120091 ) found that the radio-to-IR power-law break was 
located at frequencies higher than the if-band when GX 339-4 was in the hard state. 
However, when GX 339-4 was transitioning from the soft to hard state, the break had to 
be located at frequencies lower than H-hand as t he 0/IR data f luxes were well below that 



Gandhi et al. 



(120111 ) obtained mid-IR 



expected from the extrapolated radio spectrum. 
(WISE) with quasi-simultaneous UV (Swift), optical (Faulkes Telescope South), near-IR 
(Rapid Eye Mount Telescope) and radio (ATCA) data, and found that during the hard 
state jet emission most likely dominated i n the mid-IR and t he power-law break was also 



located in this spectral region (contrary to .Coriat et al 



20091 ) 



The hard-state optical spectrum typically remains thermal in its origin and is 
most-likely a viscously-heated accretion disk, an irradiated accretion disk and/or companion 
star, or a combination of these. 

Timing studies are valuable tools in determining the nature of emission during the 
soft and hard X-ray states. Very little has been done during the soft state in GX 339-4. 
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Motch et al.l (Il985[ ) found no sign of variability in the power-density spectrum (PSD) of 
soft state optical data. QPO's ranging between 0.05 - 1 Hz have been obser ved in the 



PSD of hard state data (JMotch et al. 



1985 : 



Gandhi 



2009 : 



Gandhi et al. 



20101 ) that may be 



attributed to oscillations within the accretion disk. Cross-correlation and auto-correlation 
functions obtained from hard state data show that optical emission l ags X-ray emission b y 



150ms (iGandhi et al. 



20081 ) while the IR lags X-rays by ~ 100ms fiCasella et al.l 120101 ). 



These timescales are much too short t o be explainec 



disk or companion star (~ 25 seconds, 



Gandhi et al. 



by re processing of X-rays by the outer 



20081 ) and suggest the dominant origin 



of hard state 0/IR emission is much closer to the black hole, possibly the base of the jet 
and/or corona. 



0/IR data taken with the ANDICAM camera on SMARTS has proven to be a va 



resource for multiwavelength studies of GX 339-4 (JHoman et al. 



2005 



Tomsick et al. 



uab 



Coriat et al. 



2009 



Buxton fc Bailvn 



Orosz et al. 



2009 



CadoUe Bel et al. 



2004 



2011 



Kalemci et al 



Corbet et al. 



2010 



20081 : 



, Dincer et al., submit t ed) and other XRB s 



2005 



Russell et al. 



Rawls et al. 



20111). 



2006 : 



Migliari et al. 



2007 : 



Using ANDICAM data of GX 339-4 



Homan et al. 



(120051 ) found strong correlations 



between X-ray and 0/IR fluxes that differed between the hard and soft states, indicating 
that the dominant flux source differed in each case. The power-law slope of the near-IR/X- 
ray correlation during the hard state was very similar to that observed for the radio/X-ray 
correlation, suggesting that the majority of near-IR emission was from a jet, although other 
components most li kely contributed to the overall 0/IR flux as well. During the soft state. 



Homan et al. 



( I2OO5I ) showed a 15-20 day lag in variations between the X-ray and if -band 
fluxes, with the if -band leading the X-ray. This time difference is consistent with a viscous 
time scale and strongly suggests that the 0/IR emission during the soft state originates 
from an accretion disk. However, since the ratio of X-ray to 0/IR flux during the soft state 
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was only a factor of 10 higher than during the hard state, it was thought the disk was either 
moderately heated or not heated at all. 



Coriat et all (120091 ) also used ANDICAM 0/IR data of GX 339-4 to discover tight 
correlations between the X-ray and 0/IR fluxes with branches that corresponded to the soft, 
intermediate, and hard X-ray states. Comparison between the IR/X-ray and optical/X-ray 
correlations suggested that the near-IR fluxes are dominated by non-thermal emission while 
the optical is dominated by thermal emission in both hard and soft X-ray states. 



Homan et al. 



(120051) 



Coriat et al. 



(120091 ) ■ and other studies have used only part of 



our 0/IR dataset of GX 339-4. In this paper, we present the entire 0/IR light curve of 
GX 339-4 from 2002-2010, inclusive, with the intent on making it public and available to 
the research community. We have also analyzed the 0/IR fluxes and SED to: 1) identify 



correlations between the near-IR and opt ical f 



such as that seen by 



Homan et al 



uxes that ar e depe ndent upon X-ray state. 



Coriat et al 



( 120091 ): 2) place limits on the 



(120051) and 

SED power-law break in the near-IR during the hard state; 3) investigate the nature of the 
flux sources during both the soft and hard X-ray states; and 4) perform a timing analysis 
of the 0/IR data to identify any characteristic timescales that may be present. In Section 
El we describe our data reduction and photometric calibration methods, and present the 
0/IR light curve of GX 339-4 and an X-ray hardness-intensity diagram (HID). In Section 
131 we distinguish two different patterns of 0/IR behavior that are associated with the 
soft and hard X-ray states. In Section HI we study the 0/IR SEDs and compare them 
to phenomenological models representing thermal and non-thermal emission. In Section 
m we investigate the variabilility timescales in the 0/IR data during both soft and hard 
X-ray states. We discuss our results in the context of other studies of GX 339-4, and give 
conclusions, in Section El 



2. Observations and Data Analysis 
2.1. Optical and Near-Infrared 



0/IR observations of GX 339-4 wer e taken with the ANDI CAMjJ camera ( iDepoy et al. 



20031 ) on the SMARTS 1.3m telescope (ISubasavage et al 



^, 



2010|). ANDICAM is a dual- 



channel imager that contains a movable internal mirror in the IR channel that allows 
dithered images to be taken while a single optical exposure is undertaken. We obtained 
pairs of images in V^ + J and I + H filters on a daily or near-daily basis whenever GX 339-4 
was available in the night sky. The optical images were single exposures of 300 seconds. 
The IR images contained 8 dithered images of 30 seconds each. The optical data were bias- 
and overscan-subtracted, and flat-fielded using the CCDPROC task in IRAF. Infrared data 
were reduced using in-house IRAF scripts that fiat-fielded, sky-subtracted, and combined 
dithered images. 

2.1.1. Magnitude Determination and Calibration 

Uncalibrated 0/IR magnitudes were determined using DAOPHOT in IRAF, for which 
the full-width-half-maximum, sky background and point-spread function parameters were 
determined separately for every image. 



Optical primary standard stars (JLandolt 



19921 ) were used to calibrate the optical 



magnitudes of three comparison stars in the field of GX 339-4. For 210 photometric 
nights we calculated the calibrated magnitude of each comparison star using photometric 
zeropoints that were measured from the SMARTS observations of the primary standard 
stars for each filter, correcting for atmospheric extinction derived from the standard star 



^http://www. ctio.noao.edu/noao/content/andicam 
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observations. The average magnitude of these three comparison stars was then used 
as a basis for differential photometry with respect to GX 339-4 for all observations. 
he near-infrared p hotometry was calibrated using near-infrared primary standard data 



( jPersson et al. 



19981 ) and two comparison stars were used as secondary standards. Only 21 
nights of photometric data were used to measure the calibrated magnitudes of the secondary 
standards as we had limited data on the near-infrared photometric zeropoints. 

We show the comparison stars in the optical and near-infrared finding charts shown in 
Figures [T] and El Our calibrated magnitudes of these comparison stars are given in Tables 
[U and [2J We also present the 0/IR magnitudes of GX 339-4 in Table [31 the full version of 
which can be accessed online. 

2.1.2. Color correction 

The {y — I) color of our secondary standards are given in Table (H and they range 
from 1.37 - 2.28. For GX 339-4, the majority of {V — I) data vary over the X-ray states 
between 0.9 and 2.0. The secondary standard colors overlap most of the {V — /) color 
range exhibited by GX 339-4. We calculated the mean color-term coefficient for the J-band 
for 1052 photometric nights between August 2003 and December 2010 to be 0.062 ± 0.048 
(l-cr). Therefore, using the average color-term correction for the /-band, the color-term 
correction for the secondary stars ranges between 0.085 and 0.141 mag. 

Color correction in the near-infrared is minimal and, as such, is not required. 
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2.1.3. Photometric Errors 

Errors on the V- and /-band differential photometry were obtained by calculating the 
1-cr error of 93 stars from 40 nights of data obtained during 2006. Results that deviated 
more than ± 3-cr from the mean were rejected, and the mean and a were recalculated 
until no more rejections were made. The l-cr error for each star was plotted against their 
calibrated magnitude and a 4th-order polynomial was fitted to the resulting distribution. 
Again we rejected any data that deviated more than ± 3-a from the polynomial and 
recalculated the polynomial fit until no more rejections were made. The same procedure 
was performed for the IR data, using 44 stars for the J-band and 29 stars for the iJ-band, 
from 40 nights in 2006. We used fewer stars than the optical bands because the field of view 
is smaller and there are fewer stars to choose from. 

Figure [3] shows 1-cr errors versus magnitude for the 0/IR and the corresponding 
polynomial fits. Once the polynomial reached a minimum, we extended that constant 
numerical value to brighter magnitudes. Thus the optical photometric errors are 0.01 mag 
for V < 17.50 and / < 16.50. Infrared photometric errors are 0.02 mag for J < 14.10 and 
H < 14.10. We refer the reader to Figure [3] to ascertain errors on fainter magnitudes. 

We also measured the error in calibrating the secondary star magnitudes by calculating 
the standard error of the mean over the number of photometric nights mentioned above. 
The photometric calibration errors are as follows: Verr = 0.003 mag, lerr = 0.004 mag, 
Jerr = 0.02 mag, Herr = 0.02 mag. These values do not account for systematic errors 
associated with such effects as the difference in effective filter responses between SMARTS 
and the standard system. Such systematics are likely to contribute a few hundredths of a 
magnitude of calibration error. Therefore we take the optical calibration errors to be ~ 0.02 
magnitudes and the IR calibration errors to be ^ 0.03 magnitudes. 

We note that in the majority of images the two nearby stars to GX 339-4 (denoted 
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as 'A' and 'B' in 



Shahbaz. Fender, fc Charles 



2OOII ) contribute negligibly to the overall 



flux of GX 339-4 during its faintest optical state as DAOPHOT is able to extract the PSF 
proflles of both GX 339-4 and stars 'A' and 'B' (the latter two are combined into one PSF 
profile). However, there are cases where the seeing or the image quality are degraded such 
that DAOPHOT was not able to separate the PSF profile of GX 339-4 from the two nearby 
stars. For completion, we have included this data in Table E] but have denoted it with an 
asterix. 

2.2. X-ray Data 

The X-ray data from proportional counter array (PGA) onboard the Rossi X-ray 
Timing Explorer (RXTE) were reduced using HEASOFT V6.7. After the standard data 
screening, the background corrected count rates were extracted from proportional counter 
unit 2 in 3 energy bands corresponding to 2.8-5.2 keV, 5.2-9.8 keV, 9.8-25 keV. 

2.3. O/IR Light Curve 



In Figure H] we present the O/IR lightcurve of GX 339-4 during 2002-2010, inclusive. 
Gaps in the data occur when the source is behind the Sun. 

There are a number of interesting features of this light curve. GX 339-4 clearly does not 
follow the pattern of long periods of quiescence punctuated by outb ursts with a fast-ris e, 



expo nential-decay morphology as is often seen in dwarf novae (e.g. 



19981 ) and neutron- star X-ray binaries, such as A ql X-1 (JMaitra &: Bailyn 



black hole binaries (JGhen. Shrader. &: Livio 



Gann 



zzo fc Mattei 



20081 ) and many 



19971 ). Rather, the source is continually active 



with strong variability, featuring rapid rises with a large range of decay times. 
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In many cases, the changes seen in the optical bands have smaller amplitudes than 
those that occur at the same time in the near-infrared bands. Obvious examples occur in 
2003 and 2005 during the rebrightening phases, discussed in Section I2.3.1I During these 
events an increase in flux is barely noticeable in the V- or /-bands, but is more obvious 
in the J- and if -bands. There are also examples of sudden drops in IR flux (e.g. mid 
2002, m id 2004 and early 2010 ) which are accompanied by more modest changes in the 



optical ( iBailyn fc Ferrara 



2004 ). These drops occur when the source transitions from the 



hard state into the soft state. This behavior is contrary to what is expected if the 0/IR 
flux originates solely frorn a thermal source, such as an accretion disk . Disk instabilities 



Cannizzo. Chen, fc Livio 



1995 



Dubus. Hameury. fc Lasota 



;asota 



200ll ) and irradiation models (JDubus et al. 



1999 



20011 ) generally attribute variability to changes in temperature 
of the accretion disk. In such cases, high luminosity should be correlated with blue color. 
However, this basic relationship clearly does not hold for GX 339-4, which suggests that a 
non-thermal emission component is often dominant. It is, therefore, desirable to study the 
0/IR SED to ascertain what flux sources are present and how they vary and contribute 
to the overall flux during these changes in the light curve. We present our 0/IR SEDs in 
Sections lU and SSI 



2.3.1. Rebrightening Events 

During 2003 and 2005, GX 339-4 displayed a slowly fading decay interrupted by a 
rebrightening which subsequently declined back, but not to the extrapolation of the previous 
decay. In Figure [5] we show these rebrightening events along with our best-flt straight lines 
to the data immediately before the rebrightening. The slopes of our flts are given in Table 
m Within errors, the 0/IR decayed at the same rate between 2003 and 2005 with the 
exception of the H band for which the decay rate had slightly increased. We extended the 
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best-fit line in Figure |5] to highlight that the decay rate after the rebrightening event is 
not a continuation of that before, but rather the 0/IR flux continues to fade to a fainter, 
steady magnitude. 

The duration of the rebrightening event in 2003 (~ 100 days) is longer than in 2005 (~ 
80 days), and the magnitude reached after the rebrightening event ended is slightly brighter 
in 2005 than in 2003 by ~ 0.5 mags. This could mean that more material was accreted by 
the black hole and/or ejected via coUimated jets in 2003 than in 2005. 

In 2007, GX 339-4 once again transitioned from the soft state to the hard state and 
experienced a rebrightening in 0/IR fluxes. However, the behavior after the rebrightening 
was markedly different to that in 2003 and 2005. GX 339-4 did not fade back to a quiescent 
magnitude but remained in the hard state during 2008 until the end of 2009, fluctuating 
significantly in 0/IR flux during that time. It is unclear why this should be the case. 

With observations of future outbursts in GX 339-4 we will explore further the 
similarities and differences between these rebrightening events and how they relate to the 
overall outburst properties. In addition, the decomposition of disk/jet flux during these 
events will be pursued in a future paper (also see Dinger et al., submitted). 

2.3.2. Hardness Intensity Diagram 



A useful to ol in discerning w'hen B HXRBs are in various states is the hardness-intensity 



diagram (HID) (JMiyamoto et al. 



19951 ) where X-ray hardness is plotted against the count 
rate. Ne arly all black- hole transients trace a 'q-shape d' path counterclockw ise in this 



diagram (JBelloni 



2005|) during it's 



20101 ). One exception is 4U 1630-47 (JTomsick et al. 
2002-2004 outburst (although it did exhibit the q-shaped pattern in 1998). In Figure |6] we 
show the HID of GX 339-4 for 2002 through 2010. HIDs of the 2002/03, 2004/05, and 2007 
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out bursts h a ve alr eady been published in 
and 



Belloni et al. 



(120051 ) 



Homan fc Bellonil (120051 ) 



Bellonil (120 lOf ). The 2008/09 data are an extension of the outburst started in 2007. 
The 2006 and 2007-09 outbursts are interesting for the following reasons. In 2006, GX 
339-4 experienced an outburst relatively short compared to the others and that never left 
the hard state. In the HID we see this outburst only resides in the lower vertical branch of 
the 'q' and never resides in the loop. However, we see in Figure S] that the 0/IR flux did 
increase significantly but did not peak at the same fluxes as in 2002, 2004, 2007 and 2010. 
As already noted, the 2007-2009 outburst remained in the hard state for a very long time 
after transitioning from the soft state, and in the 0/IR we see it fluctuating signficantly 
during this time. The HID for this outburst does not show any obvious differences to those 
of other outbursts. The only difference being the amount of ti me it has spen t on the lower 



horizontal branch, also known as the hard-intermediate state 
outburst looks to be typical of other outbursts in the HID. 



Bellonil tom . The 2010 



3. O/IR States in GX 339-4 



3.1. O/IR Correlations 



Coriat et al.l ( 120091 ) described correlations between the X-ray and V- and if -band 
fluxes in GX 339-4 in which a clear delineation of data is seen between the hard, soft, and 
intermediate X-ray states. Here we investigate correlations within the O/IR and show how 
this relates to the X-ray behavior. 

In Figure [7] we plot the K-band magnitudes against J-, J-, and H-hand magnitudes. 
In each plot there are clearly two separate branches, and the separation between the two 
branches increase from the I to H band. The slope of the upper branch seems to change 
near the brighter end. To confirm this, we performed both single and broken power-law fits 
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to each branch using the Levenberg-Marquardt algorithm. The best-fit parameter values 
are summarized in Table O We found the lower branch was best described by a single 
power-law with a slope less than one, that is, GX 339-4 became bluer as the system became 
brighter. The slope of the lower branch decreased from the I to H band. The upper branch, 
on the other hand, was best fit by a broken power-law. From the faintest magnitudes to the 
break, the slope is greater than one indicating that the system gets redder with increasing 
brightness. This slope increases from the I to H bands. Beyond the break, the power-law 
slope is less than one, as in the case of the lower branch but with a higher slope in the J and 
H bands than that observed for the lower branch. The location of the break is statistically 
indistinguishable (within 3-cr) between the V/J and V/H correlations, but is located at a 
slightly fainter magnitude in the V/I correlation. 

The color vs. brightness trends suggest that, from the faintest optical state up to the 
break, the upper branches represent a state in which the 0/IR is dominated by changes 
in non-thermal emission while at magnitudes brighter than the break, and in the lower 
branch, suggest changes in a thermal component is more dominant. Since it is regarded 
that non-thermal emission is from collimated jets, the 0/IR flux during the hard state 
must contain components associated directly with the inner accretion flow in addition 
to whatever flux is produced directly by the outer accretion disk and the (much fainter) 



companion star. 



Maitra fc Bailvnl toOSJ ) and 



Russell et al 



( 120111 ) presented analysis of the 



0/IR color-magnitude diagram for Aql X-1 and XTE J1550-564, respectively, and showed 
that the data could be well explained by a heated blackbody (representing the accretion 
disk) plus non-thermal jet emission during the hard state. These results support the above 
interpretation of an accretion disk dominating the soft state, and jet emission dominating 
the hard state. In a future paper we will present similar analysis of color magnitude 
diagrams for GX 339-4. 
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Association between th e correlation slopes and the jet/disk emission could already be 



inferred from the results of 



Coriat et al 



(I2OO9I ). who explored the relationships between 



0/IR flux and X-ray flux and ascertained that the upper branch corresponds to the "hard" 
X-ray state, and the lower branch to the "soft" state. We conflrm the association between 
the 0/IR upper and lower branches and the X-ray hard and soft state as follows. First, we 
deflned the "soft" and "hard" X-ray states using the RXTE PGA hardness ratio (HR) data 
(where HR is the ratio of fluxes in the bands (9.8-25 keV)/(2.8-5.2 keV)). The entire GX 
339-4 RXTE dataset was analyzed to obtain count rates and hardness ratios. The beginning 
and end of outbursts in the hard state were determined using the light curve and setting 
a threshold of 10 cts/s/PGU. Then we created spectra and power density spectra in the 
3-25 keV band for all observations with good time intervals exceeding 1 ks while the count 
rate increased. We stopped the analysis when signiflcant changes were observed either in 
the X-ray spectra or the power spectra indicating a transition to an intermediate state. 
For all these observations, the photon index was less than 1.75, and the power spectrum 
was dominated by broad Lorentzians with mean rms amplitudes greater than 25%. The 
hardness ratio threshold including all observations was 0.6, however, for some outbursts the 
threshold can go up to 0.75. When categorizing the 0/IR into soft and hard states, we 
defined the soft state to be when HR < 0.6, and the hard state when HR > 0.6. Second, we 
considered data only for days when simultaneous V- and if-band data were taken, as this 
pair of data show the separate branches most clearly. Using these criteria, we plotted in 
Figures |8] and |9] red data points when the HR < 0.6 (soft state), black points when HR > 
0.6 (hard state), and blue points when the X-ray state could not be determined due to low 
X-ray counts. We show all data as green points in Figures [8] and [9] for completeness. The 
analysis clearly shows that the lower (red) branch corresponds to when GX 339-4 is in a 
soft X-ray state, while the upper (black) branch corresponds to the hard X-ray state. This 
confirms that we are able to use 0/IR data alone to determine when GX 339-4 is in either 
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the hard or soft X-ray state. 



3.2. State Transitions 



CoriatetaL 



(120091) discussed how GX 339-4 moved along the IR/X-ray correlation 
branches over the course of four separate outbursts (2002/2003, 2004/2005, 2006 and 
2007). In the following sections we trace the paths of the state transitions in the 0/IR 
correlations and show the hard-to-soft transition in 2010. We give a timescale over which 
the hard-to-soft transition occurs in 2002 and 2010. 



3.2.1. Hard-to-Soft Transition 



In Figure HOI (left) we show the two best examples of transitions from the hard to soft 
X-ray states during 2002 (red line) and 2010 (green line). These are the two outbursts 
during when we had the best sampling over the hard-to-soft transition (i.e. at least one 
observation per day with no gaps between days over the transition period). When GX 
339-4 begins a new outburst (that is from the faintest optical state) it is in the hard state 
at the lower end of the hard-state branch and rises up this branch to the brightest end. 
This corresponds t o the far right vertical branch in the HID (Figure ED, also known as 



the low-hard state 



Bellonil (J2010l ). Some time after GX 339-4 reaches its peak brightness 



it begins to transition to the soft state (the upper horizontal branch in the HID = 
hard- intermediate state). During this p h ase G X 339-4 rapidly fades, more in the IR than 



the optical. We note that JMotch et al 
brightness over the transition 
during 2002 was also noted by 



( Il985l ) was the first to see a change in optical 



Tom hard to s oft st ates. This sudden drop in 0/IR flux 



Homan et al. 



(120051). In both 2002 and 2010 the path taken 



by GX 339-4 over this transition is remarkably similar. Once it has left the hard-state 



-18- 

branch it takes around 4 days to reach the soft state branch. Both cases merge with the 
soft branch at arou nd the same location: V ~ 16.5 mag. From our data, and as stated in 



Coriat et al. 



(|2009[ ). we can only say that the rise of GX 339-4 follows the same path as 
the decay phase (see Section I3.2.2p from a V^-band magnitude of around 18.2. We do not 
have any data to state definitively that the rise and decay paths are the same below this 
magnitude. 

3.2.2. Soft-to-Hard Transition 

When GX 339-4 is in the soft X-ray state (the left vertical branch in the HID = 
high-soft state) it is already at relatively bright 0/IR fluxes. In Figure ^U\ (right) we show 
three examples (2003, 2005, 2007) of a soft-to-hard transition (corresponding to the bottom 
horizontal branch in the HID). During 2003 and 2005 the transition from soft to hard X-ray 
states followed the same path: GX 339-4 moves down the soft state branch until it reaches 
the intersection between the soft and hard state branches. It then moves up the hard state 
branch to a peak magnitude (slightly brighter in 2005 than 2003 by around 0.2 mag in 
V) and then swiftly down the hard state branch to the faintest end. After transitioning 
into the soft state in 2007 GX 339-4 did not move down to the faintest end of the hard 
state branch until 2009. Throughout this period GX 339-4 moved up and down the middle 
portion of the hard state branch. Toward the end of 2009 it finally moved down the hard 
state branch to the faintest end. 
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4. Spectral Energy Distributions 



4.1. Determination of O/IR Flux 



To obtain fluxes appropriate for comparing with observations in other wavelength 
regimes requires two steps. First, the magnitudes must be dereddened so that they are not 
affected by absorption along the line of sight; second, the dereddened magnitudes must be 
converted to flux units. 



19981 ) that was converted to Ay = 3.7 ± 0.3 (jCardelli. Clayton, fc Mathii 



The 0/lR magnitudes were dereddened u sing EfB-V) = 1.2 ib 0.1 (IZdziarski et al. 

id 



19891). We 



calculated the extinction coefficient in the other ba nds using A(X)/A (V) = a{x) + b{x)/Rv 



and the relationships for a(x) and b(x) as given in 



and 



Cardelli. Clayton. &: Mathis 



km) 



O'Donnelll fll994h for V and / bands. 



for J and if bands. In calculating a(x) a nd 



b(x), we used the effective wa velength for each band (JFrogel et al. 



1978 



Elias et al. 



1982 



Bessell. Castelh. fc Plez 



1998h : V = 545 nm, I = 798 nm, J = 1250 nm, H = 1650 nm. The 



errors on the dereddened magnitudes are dominated by the error in Ay, and are as follows: 
Verr = 0.3 mag, lerr = 0.1 mag, Jerr = 0.1 mag, Herr = 0.1 mag. Adding the photometric 
and interstellar redenning errors in quadrature give the following total errors on dereddened 
magnitudes: Verr = 0.30 mag, lerr = 0.10 mag, Jerr = 0.10 mag. Hem = 0.10 mag. 



We note that 



Roman et al. 



(120051) found an independent value of E(B-V) using 



iVf/ = 5 ± 1 X 10^^ atoms cm'^ fiKong et al.ll2000h . convert ing to EfB-V) = 0.94 ± 0.19 



Predehl &: Schmitt 



19891). 



Zdziarski et al 



1995h. in turn giving Ay = 2.9 ± 0.6 fICardeUi. Clavton. fc Mathis 



(119981 ) obtained their value of E(B-V) from a detailed study of 
extinction of ~450 stars within ±5° of GX 339-4. In either case, E(B-V) agree within the 
errors. 



To convert the dereddened magnitudes to flux (in units of Jy) we used the following 
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zeropoint fluxes (where the zeropoint flu x of a given filter is that 



magnitud e): Vo = 3636 Jy, In = 2416 Jy ( 



980 Jy flFrogel et al. 



1978 



Elias et al 



Bessell. Castelli. fc Plez 



corre sponding to zero 



19981), Jo = 1670 Jy, Hq 



19821 ). The dereddened fluxes of GX 339-4 for each 



waveband are given in Table O 

We have checked our flux calibration method by obtaining observations of optical 
spectrophotometric standar ds (LTT 3218 , LTT 4 816, EG 131, EG 27 4) an d compa r ing o ur 



fluxes to those observed by 



Hamuy et al. 



(119921 ) 



Hamuy et al 



fll994J ). and 



Besselll fll999[ ). 



Observations were obtained in B, V, R, and /-bands and our magnitudes were converted 
into fluxes using the method outlined in Section 12.1.11 We compared the photometric fluxes 
to the flux calibrated spectrum over the relevant bandpass, and found no deviations larger 
than 0.05 dex, or 10% of a given flux. We show two examples (LTT 4816 and EG 274) in 
Figure HH 

4.2. Soft state SED 



As stated in Section 13. H our flts to the slope of the soft-state branch indicate that the 
flux increases more quickly at bluer wavelengths. In other words, the source gets bluer as 
it gets brighter. Figure [12] shows two example SEDs obtained from the soft-state branch 
data for V = 16.0 and 17.0 mag. The J, H and /-band data were extrapolated from the 
power-law flt to the correlation data shown in Figure [71 In both cases, the 0/IR data 
lie on the same spectral slope. We determined the slope by fltting a 1-deg polynomial in 
log-log space using the Levenberg-Marquardt method and derived the following results: 
a = 1.6 ± 0.1 for V = VJ mag; a = 1.8 ± 0.1 for l^ = 16 mag. Hence, the slopes in each 
case are the same, within errors. We show these flts as dotted lines in Figure [T^l Also 
shown, as a dashed line, is the slope expected for F oc v^/^ representing a viscously-heated 
steady-state accretion disk. The 0/IR SEDs during the soft state have spectral slopes closer 
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to that expected for the RJ tail [a = 2). An RJ tail may be present from either the outer 
regions of a viscous-heated disk or an irradiated disk. Since we do not observe a spectral 
turnover we cannot say definitively what the temperature of the thermal source is nor what 
the heating mechanism is. In any case, the trend seen during the soft state in GX 339-4 
with higher fluxes is consistent with a thermal source that is heated to higher temperatures. 

4.3. Hard state SED 

In Section 13.11 we found that the hard state branch was best described by a broken 
power-law where at fainter magnitudes GX 339-4 gets redder at higher fluxes while beyond 
the break the opposite occurs. For the former case, this requires either a non-thermal source 
or a thermal source that dramatically increases in size as it decreases in temperature. For 
the latter case, GX 339-4 is getting hotter with higher fluxes. 

Figure [13] shows the evolution of the hard-state 0/IR SED from the faintest V-band 
magnitude {V = 19.5) to the brightest observed on the hard-state branch {V = 15). At 
the faintest magnitude, the 0/IR data lie on the same spectral slope {a = 1.3 ±0.1). This 
is shown as the solid line in Figure [131 This slope is less than the slope observed during 
the soft state. The spectral slope expected for a viscous-heated accretion disk {a ^ 1/3) is 
also shown as the dashed line. As for the soft state, we see that the optical bands have a 
much steeper slope but we are unable to discern the temperature or heating mechanism of 
the thermal source. Nevertheless, even when GX 339-4 is at its faintest optical state we see 
evidence of a heated thermal source, most likely the accretion disk. In Figure [13] we have 
propagated the fitted slope at the faintest magnitude to the SEDs at brighter magnitudes to 
show that, as GX 339-4 increases in brightness, the V and I bands are consistent with this 
slope, indicating that V and / are most likely domi nated by a the r mal so urce that is heated 



to higher temperatures with increasing 0/IR flux. 



Homanetal 



(12005! ) also reported a 
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stable spectral slope between the V and / bands during the hard state rise and transition 
to the soft state during the 2002/03 outburst. 

The near-IR bands, on the other hand, show a clear evolution in slope that decreases 



with increas ing flux, whi 



observed by 



Homan et al. 



e the 



flux in the near-IR bands also increases. This was also 



fl2005h during the 2002/03 outburst. 



Radio emission in GX 339-4 and other XRBs has only bee n observed d ur- 
ing the hard state and transitions from one state to another ( 1 



Corbel fc Fender 



2004 



2002 



Corbel et al. 



Fender. Homan. fc Belloni 



2003; 



Fender 



Fender. Belloni. fc Gallo 



2004 : 



2001 



Gallo et al. 



20091 ) . Presumably the radio emission is due to syn- 



chrotron emission from a jet. It has been suggested that such emission might contribute 



signiflcantly at near-IR bandpasses ( ICorbel fc Fender 



2002 



Nowak et al. 



20051 ) and explain 



the near-infrared excess we see in our hard-state SED. We searched the literature for radio 
observations that were taken on the same night as the SMARTS data during the har d state. 



20051 ) with 



We found two cases. The flrst was taken on UT 2002 April 03 (iNowak et al. 

ATCA at 1.38 GHz (4.83 ± 0.20 mJy) and 2.40 GHz (5.13 ± 0.11 mJy). On this date, we 

observed GX 339-4 with a \ ^-band magnitude of 15.73. The second case was on UT 2002 



April 18 f.Nowak et al. 



2005|) at 8640 MHz (13.4 9 ± 0.08 m, 



mJy). Note 



t hat th is data was flrst reported in 



Nowak et al 



Gallo et al 



y) and 4800 MHz (12.97 ± 0.07 



( I2OO4J ) but was re-analyzed by 



(120051 ). On this date, GX 339-4 had a l^-band magnitude of 15.14 mag. We 



show these radio data along with our 0/IR hard-state SEDs for these \^-band magnitudes 
in Figure [HI We performed power-law flts to 1) the radio, J and H band data; and 2) 
radio data only, using a deg-1 polynomial and the Levenberg-Marquardt method. These 
flts are shown in Figure [T3] as the solid and dotted lines, respectively. The best-flt slopes 
are for SED 1 (03 April): aradio/m = 0.14 ± 0.01; aradio = 0.11 ±0.09; and for SED 2 (18 
April): aradio/iR = 0.14 ± 0.01; aradio = 0.14 ± 0.01. 
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On UT 03 April 2002, the error on the radio data is sufficiently large that a simple 
extrapolation from the radio to the near-IR spectral may or may not be consistent with 
the 0/IR data. A fit to the radio and near-IR data on this date has the same spectral 
slope as that found on UT 2002 April 18. The latter has much tighter constraints thanks 
to the relatively small errors on the radio data. Indeed, the spectral slopes are the same 
with and without the near-IR data included in the fit. Interestingly, the J-band also lies 
on this extrapolation. Since the /-band lies on the same spectral slope as the l^-band 
from the faintest to brightest magnitudes during the hard-state (i.e. the V — I color is 
constant) then the I band is most likely the location of a break between the non-thermal 
and thermal SEDs. We do not observe a break in the non-thermal spectrum, however, we 
can say that the non-thermal bre ak location m u st exi st at frequencies greater than or equa l 



to the J-band, in agreement with 



Coriat et al. 



(120091 ) but contrary to 



Gandhi et al. 



(1201lh 



(see Section [6] for further discussion of this point). 

At magnitudes brighter than the break in the hard-state branch, GX 339-4 gets bluer 
at higher fiuxes. The SED shown in Figure [13] shows that the non-thermal component is 
present, even at the brightest magnitudes (i.e. beyond the break). However, the near-IR 
fiux at V = 15.0 mag relative to the V-I spectral slope is less than at V = 16.5 (just before 
the break). In other words, as the V-I fiuxes increase beyond this break, the near-IR fiuxes 
increase, but less so than at magnitudes below the break. Clearly either a physical change 



has occurred in the 
339-4 gets brighter. 



et and /or th e 



Coriat et al. 



jet o utput may be reaching a saturation point as GX 



( 120091 ) suggest that the break in the jet spectrum SED 
moves to higher frequencies at higher fiuxes. Below the break in the hard-state branch we 
may be observing optically thin jet emission while above the break at higher fiuxes we may 
be observing optically thick jet emission. We leave further investigation of this point to a 
future paper. In addition, changes in the thermal emission become more dominant beyond 



the break. Since the V — I slope is the same this simply translates to the thermal source 
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(i.e. disk) getting hotter with increasing flux. 

In summary, simultaneous radio and near-IR data can be explained by a non-thermal 
component, and this component dominates at all but the faintest magnitudes. The optical 
spectral region remains dominated by a thermal flux source, most likely a heated accretion 
disk. Therefore, the 0/IR SED needs to be explained by a combination of thermal and 
power-law models. We are unable to perform such a flt since the number of parameters 
involved would exceed the number of constraints (data) available. This highlights the need 
to obtain simultaneous data at other wavelengths, such as B, UV, and mid-IR in addition 
to the 0/IR, radio and X-ray, to further sample the shape of the SED near the 0/IR and 
to best constrain model flts. In addition, obtaining simultaneous radio data with small 
errors and/or as many bandpasses as possible would strengthen model flts and place more 
confldence in associating near-IR flux with jet emission. 

5. Timing 

We have calculated the power spectral density (PSD) of the variability of GX 339-4 in 
the hard and soft states. The PSD corresponds to the power in the variability of emission as 
a function of timescale. The power in the units of rms^ Hz~^ is given by: P(z/)=Az/°, where 
A is a normalization constant, u is the frequency (inverse of timescale) in Hz and a is the 
power law index. In a log[P(z/)] vs. log[z/] plot, a is deflned as the slope of the power-law 
power spectrum. The timescales covered in the power spectra range from ~8 months (the 
longest well-sampled intervals during which any state transition did not take place) to 1-2 
days (the most frequent sampling carried out for a signiflcant interval). Hence, the power 
spectra range from 10^^'^ Hz to 10^^'^ Hz. 

We divided the V- and //-band light curves into "soft" and "hard" intervals. 
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respectively, based on the hardness ratio in the contemporaneous RXTE-PCA observations, 
as outhned in Section |31 These two bands were chosen as the H-hand data shows hght-curve 
variabihty most clearly during the hard state, while the same is true for the \^-band during 
the soft state. However, PSDs were calculated for each band in both the hard and soft 
states. The time intervals chosen were V: JD = 2453224.70 - 2453478.79; and H: JD = 
2454498.88 - 2454746.55. These were the longest well-sampled intervals during which the 
X-ray hardness ratio stayed soft and hard, respectively. Over these time intervals, the 
V-band magnitude ranged between 15.9 and 17.6 mags, while the if -band varied between 
12.7 and 15.3 mags. 

The raw PSD calculated from a light curve combines two aspects of the dataset: (1) the 
intrinsic variation of the object and (2) the effects of the temporal sampling pattern of the 
observations. In order to remove the latter, we applied a M onte-Carlo type algorithm based 



on the "Power Spectrum Response Method" (PSRESP) of 



Uttley. McHardy. fc Papadakis 



( 120021 ) to determine the intrinsic PSD (and its associated un c ertain ties) of the light curves. 



Our realization of PSRESP is described in iChatterjee et al. 



(120081). PSRESP gives both 



the best-fit PSD model and a "success fraction" -Fgucc (fraction of simulated light curves 
that successfully represent the observed light curve) that indicates the goodness of fit of the 
model. We bin the data in 1 to 2-day time intervals, averaging all data points within each 
bin to calculate the flux. We filled empty bins through linear interpolation of the adjacent 
bins in order to avoid gaps that would distort the PSD. We accounted for the effects of the 
binning and interpolation by inserting in each of the simulated light curves the same gaps 
as occur in the actual data and performing the same binning and interpolation procedures. 

In Figure [15] we show the PSDs of GX 339-4 and the best-fit model in H (hard state) 
and V bands (soft state). The PSDs show red noise behavior, i.e., there is higher amplitude 
variability on longer than on shorter timescales. Based on the model with the highest 
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success fraction, the H-hand PSD, and therefore the hard-state PSD, is best fit with a 
simple power law of slope — LSI^q'^, for which the success fraction is 0.56. That for the 
V-band (soft state) PSD is — 2.0J;q'5 with success fraction 0.89. During the fitting, we 
varied the slope from —1.0 to —3.0 in steps of 0.1. The uncertainties of the slopes represent 
the range of slopes beyond which the success fraction goes below half of the maximum 
value. The rejection confidence, equal to one minus the success fraction, is much less than 
0.9 in both cases. This implies that a simple power-law model provides an acceptable fit 
to the PSD at both wavebands. The PSD for the V band during the hard state, and H 
band during the soft state, were also well-fitted by a simple power law model with slopes 
consistent with those described here within uncertainties. A summary of our PSD fits are 
given in Table [6l We see no break in the power spectra, nor evidence of any QPO's, at few 
days to few months timescales. 

The standard deviations of the light curves are 7.8 mJy (l^-band) and 3.6 mJy {H- 
band). The mean fluxes are 27.6 mJy (\^-band) and 6.2 mJy (if -band). The uncertainties 
are negligible compared to the variation. To meaningfully compare the fluctuation of two 
light curves, we need to i) normalize the standard deviation by the mean of the light 
curves ii) account for the contribution of the observational uncertainties, as opposed to the 
intrinsic fluctuations, to the observed variability of the light curves. Here, we calculated the 
"excess variance" normalized by the square of the mean flux, as a measure of the variability 
in the V- and if -band data. It is deflned as 



K 



b (T„,., 



var \ / 9 



where S is the variance, ( Xprr is the observatio i ial uncertainty, and x is the mean of the 



data (JNandra et al. 



1997 : 



Edelson et al. 



2002 



Vaughan et al. 



20031). We flnd that the 



values of i^var of the V- and ii-band variability for the entire data set are 0.28 and 0.58, 
respectively, i.e, of the same order of magnitude with the ii-band having a higher value. 
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This is consistent with the shghtly larger amphtude of power in the H- compared to the 
l^-band as seen in Figure [T3 

The timescales covered in Casella et al. (2010) and Gandhi et ah (2010) are from 
few minutes to 0.1 sec (0.003 Hz to 10 Hz). Therefore, there is a large gap between the 
frequency range covered in those papers and that presented here. The slope at the lowest 
frequency end of the power spectrum shown in those work is significantly flatter than the 
slope determined here for either of the wave bands. But the data used in this paper is 
insufficient to determine at which frequency the flattening takes place. 

Well-sampled light curves covering minutes to days timescales can be used to fill-up 
the gap in the frequency range and check the broadband nature of the power spectrum from 
months to sub-second timescales. This may indicate whether the slope changes gradually 
over a range of timescales, or it changes abruptly and in the latter case, the numerical value 
of that characteristic timescale. This will be carried out in a future paper. 

6. Discussion and Conclusions 

We have presented a comprehensive optical and near-infrared light curve of GX 339-4 
between 2002 and 2010, inclusive. The light curve shows that GX 339-4 is continuously 
in outburst (that is, the companion star never dominates the 0/IR fiux) but experiences 
significant changes in luminosity and spectral state. 

When GX 339-4 transitions from the hard state to the soft state changes in the near-IR 
fiuxes are much greater than in the optical fiuxes. This can be understood if the near-IR 
fiux is dominated by non-thermal emission from, say, a jet. The change in near-IR fiux is 
significantly greater during the hard-to-soft transition than the soft-to-hard transition (i.e. 
rebrightening event) suggesting that non-thermal emission is more dominant during the 
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hard state before the hard-to-soft transition than it is after the soft-to-hard transition. 

We have seen three examples of rebrightening events during 2003, 2005, and 2007. 
However, the 2007 event was markedly different to the others in that GX 339-4 did not 
return to quiescence after the peak was reached, and remained in the hard state during 
2008 and 2009 during which it fluctuated significantly in 0/IR fluxes. Rebrightening eve nts 



have been observed in the 0/IR during th e outbursts of 4U1543-47 (JBuxton &: Bailyn 



2004) and XT 



( JRussell et al. 



J1550-564 flJain et al 



20011 ). and in optical/X-ray fluxes in XTE J1752-223 



20121 ). The shape of the rebrightening events in all three systems are 
remarkably similar in that they exhibit a fast-rise-exponential-decay profile. The increase 
in 0/IR occurs soon after the source transitions from the soft to the hard X-ray state (via 



the inter mediate state in 4U 1543-47, 



emission (ICorbel et al. 



2001 



Kalemci et al. 



Jain et al. 



2001 



2005) 



Kalemci et al 



ong w ith a resurgence in radio 



20051 ). The source then decays 



back to a magnitude where it remains steady until the next outburst (except for the 2007 
event in GX 339-4). Since: 1) these events occur after transition into the hard state; 2) 
the flux becomes redder with increasing flux; and 3) radio emission is detected near the 
peak of the 0/IR maximum, it seems likely that rebrightening events are associated with 
the reappearance of collimated jets. The rebrightening events observed in GX 339-4 have 
differences between their durations and the final magnitudes reached after the peak has 
faded. Only one rebrightening event has been observed in 4U 1543-47, XTE J1550-564, and 
XTE J1752-223 so we cannot comment on any differences that may occur from one outburst 
to the next in these systems. However, the duration of the rebrightening events in other 
BHXRBs are significantly different to that of GX 339-4: ~ 60 days in XTE J1550-564, ~ 
35 days in 4U 1543-47, and ~ 40 days in XTE J1752-223. These differences may provide 
clues as to how much material has been either accreted by the black hole or ejected via the 
jets. Research is now starting to focus on ho w much iets contribute to the overall flux of 



the system at high and low luminosities (e.g. 



Russell et al. 



2010 



Soleri fc Fender 



20111 ) and 
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decomposition of the disk/jet fluxes (e.g. Maitra et al.ll2009l : iRussell et al.ll2010l . I2OIII l2012l . 
Dinger et al., submitted). For GX 339-4, we are now able to observe multiple rebrightening 
events for one BHXRB. The differences we have observed in the rebrightening events may 
provide interesting and crucial limits to theoretical models, and we may now begin to study 
whether these differences are correlated with other properties of the outbursts. 

The HID of GX 339-4 is broadly similar for all outbursts except in 2006 when GX 
339-4 failed to travel around the HID loop. Based o n an emp i rical r elati onship between 



the hard X-ray peak flux and outburst waiting time, .Yu et al. 



(120071 ) and 



Wu et all teoioh 



suggest that, assuming that the hard X-ray flux produced is directly related to the amount 
of mass accreted from the disk, there exists a disk mass limit (as yet unknown) that must 
be exceeded in order for GX 339-4 to generate a hard X-ray outburst above ~ 0.2 Crab. In 
2006, GX 339-4 only reached ~ 0.1 Crab, below this limit. Future studies may shed light on 
what this limit is physically and quantitatively, and how it is related to the HID. GX 339-4 
experienced an outburst in 2010 that was above this limit. However , it occurr e d som ewhat 
later and had a lower peak flux (~ 0.4 Crab) than that prediced by 



Wu et al. 



torn . The 



other outbursts move around the loop counterclockwise, however there are some differences 
1) between the upper and lower horizontal branches, especially during 2004/05; 2) the 
extent to which GX 339-4 reaches the far left part of the diagram; and 3) the amount of 
time spent in each branch of the HID. We may now be close to a position that allows us to 
start studying correlations between the various properties in the HID (X-ray) and the 0/IR 
outbursts. 

We have found strong correlations within the 0/IR with two branches corresponding 
to either the hard or s oft X-ray states. T hese correlations are similar to those seen between 



0/IR and X-ray flux (Coriat et al. 



2OO9I ) with distinct branches corresponding to the 



hard and soft states. With the correlations so clearly deflned, particularly between V- 
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and if -band, it is possible to determine whether GX 339-4 resides in the hard or soft 
X-ray state from the 0/IR data alone. We have found the soft state branch is best 
described b y a single power-law while the hard-state branch is best represented by a broken 



power-law. 



Coriat et al. 



( 120091 ) also found a break in their hard state branch, but only in 
the if -band/X-ray correlation (not in the V-hand/X-iaj correlation) and at significantly 
lower ii-band fluxes than that observed in this work. Their break occurs just before the 
the hard-state branch bifurcates whereas we observe a break near the brightest end of the 
hard-state branch. In future work, we will explore the difference in break location and how 
this break may or may not relate to the properties of the jet spectrum. 

In this paper we showed the paths traced in the 0/IR correlation space when the 
hard-to-soft and soft-to-hard state transitions occurred. For the former, we showed two 
well-sampled examples (2002, 2010) and that this transition occurs very quickly (~4 
days). The paths taken in these two cases were very similar and ended the transition in 
the soft-state branch at around the same magnitude: V ~ 16.5 mag. For the latter case 
we showed three cases. In two (2003 and 2005) we saw very similar behavior where the 
transition ended at the faintest end of the hard state branch. In the third (2007), however, 
it moved up and down the hard-state branch (never to the faintest end) where it remained 
during 2008/2009. In future observations of GX 339-4 we will closely monitor the state 
transitions in the 0/IR to see whether the paths traced in the correlation space vary or 
otherwise with each outburst. In particular, we will see if the rise during the hard state 
from the faintest optical state to the brightest follows the same path as the decay phase 
at the faintest end of the hard-state branch. It will also be interesting to see whether the 
location of the soft-branch merge remains the same in future hard-to-soft transitions. 

The soft state 0/IR SED is consistent with a thermal flux source that is most likely 
a heated accretion disk. We were unable to discern whether this source is irradiated or 
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not. The spectral slope between V and / bands during the hard state is a = 1.3 ± 0.1. 
For the soft state a = 1.6 — 1.8 ± 0.1, slightly h igher than the hard-stat e optic al slope, 



and consistent with thos e found by 



BHXBs (0.5 < a < 1.5). 



Coriat et al. 



HvnesI (120051) and 



Zurita Heras et al. 



(1201 ih in other 



(120091 ) measured the spectral slope between the H- 
and l^-band data to be 1 < a < 2 during the soft X-ray state and found it to be positively 
correlated with V"-band flux. We also see the 0/IR slope increasing with increasing 
brightness during the soft state, however, they are consistent within our errors. 

The hard-state SED shows evidence of both a non-thermal and thermal component. 
A power-law spectrum extrapolated from the radio data in the hard-state SED during the 
2002 outburst is consistent with the near-infrared data but not the optical. The slopes 



derived 



or the radio /near-IR spectrum (a = 0.14 — 0.17 ± 0.1) are consistent with other 



studies (1 Corbel &: Fender 



2002 



Homan et al 



20051 ). On 18 April 2002 the errors on the 



radio data are sufficiently small to allow us to place tight constraints on the non-thermal 
power-law slope. We do not see evidence of a non-thermal power-law break at frequencies 



lower than the J band, suggesting that the brea 
equal to the J b and. This is in agree ment with 



^ is located at freq uencies greater than or 



Coriat et al. 



was observed by 



Gandhi et al. 



(120091 ) but contrary to what 



(1201 ll ). During 2010 when GX 339- 4 was once again very 



Gandhi et al. 



(I2OIII ) found the power-law 



bright during the hard state {V = 14.74 mag) 
break to be located in the mid-IR. They fit a power-law to the radio- W4 data with a slope 
of +0.29 ± 0.02, significantly higher than slopes found in this paper and others. Since we 
do not have simultaneous mid-IR data for the 2002 outburst, we cannot comment on what 
differences or similarities may exist in this spectral region. However, if one considered just 
the 2010 radio data alone, the errorbar on the 9 GHz is sufficiently large that the radio 
spectrum could lie well below the mid-IR data and line up with the near-IR instead. As 
Gandhi et al. also note, fitting single power-laws in the radio-to-near-IR spectral region 
may be an oversimplification. Using sophisticated theoretical models, such as those shown 
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in iMarkoff et al.l ( 120051 ) , iHoman et al.l ( 120051 ) , and iMaitra et al.l (120091 ) may be necessary to 
explain the connection between the radio/mid-IR/near-IR data. 



Coriat et aL 



(J2009[ ) used as yet unpubhshed radio and our SMARTS 0/IR data to 
construct an SED when GX 339-4 was transitioning from the soft state to the hard state 
(UT 2005 Apr 21, MJD 53482). The 0/IR fluxes lie well below the radio spectrum of slope 
a = 0.13± 0.02. Eleven days later (MJD 53493), the radio spectrum lines up with the 
near-IR with a slope a = 0.17 ± 0.01. At this time, GX 339-4 reached its peak in the 
IR flux while the radio flux was declining. Nine days later (MJD 53502) the radio flux 
decreased further while the 0/IR flux decreased slightly and the near-IR was consistent 
with a radio spectrum of slope a = 0.21 ± 0.02. Hence the power-law spectrum was 
significantly different during a soft-to-hard state transition than during the hard state 
during a rebrigtening phase. Note this was at fainter 0/IR mags than our SEDs in Figure 



To model the hard-state 0/IR SED would require a greater number of parameters than 
the data we have at hand can constrain. A s noted above, more s ophisticated model-fi ts of 



the hard-stat e SED have been perfo r med 



Maitra et all ( l2009h . 



Markoff et al. 



by 



Markoff et al. 



(120051 ) 



Homan et al. 



(120051), and 



( I2OO5I ) included data fror n the radio and X-ray regimes 



and u sed q uasi-simult a neous 0/IR data obtained in 1981 by 



(ligSlh and 



Pederson 



Motch. Ilovaisky. &: Chevalier 



(Il98ll ). Their concl usions were that the 0/IR could be explained 



by a combination of a multi-temperature 



Shakura &: SunyaevI (119731 ) disk and a single- 



Homan et al. 



(I2OO5I ) applied models to 



temperature blackbody representing the irradia tion of the ou t er ac cretion disk by the 
jet, with the IR dominated by fiux from a jet. 
simultaneous SMARTS 0/IR (V, /, and H bands), radio and X-ray data taken during the 
2002 outburst. Two cases were considered: a pure-jet model (that is, no contribution to 
the optical by a disk), and one where the V and / band flux originate from a multicolor 
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thermal disk and the I to H bands from a jet. In both cases, th ere were difficu 
fitting the models to the steep V-to-I band slope. We note that in 



Homan et al 



ties in 



fl2005[ ) the 



/-band was not dereddened correctly. In our paper we have corrected this and, as a result, 
we would see the J-band move up in the 0/IR SED, producing a shallower slope between 
the V and /bands that would ease any difficulties Ho man et al. had in fit ting the models. 



Maitra et al. 



(I2OO9I ) also used the same 0/IR data as iHoman et al.l (120051 ) and concluded 



that their radio/0/IR/X-ray SED could be entirely explained by a jet-only model. In light 
of the /-band dere ddening correc t ion, w e str ongly encourage a re -analysis of such model 



fitting. In addition 
(V-, / -, and //-band 



Homan et al 



Homan et al. 



2005) and 



Maitra et al. 



(120091 ) used only three points 



20051 ). whereas the 1981 dataset used by 



Markoff et al. 



(120051 ) included 9 data points from B to the near-IR. Had data at shorter wavelengths 



been available, Hom an et al. and Mait r a et a l. may have seen evidence for higher flux at 



shorter wavelengths (iCadoUe Bel et al.l (1201 ll ) see a rise in flux from optical through to UV 



wavelengths) which may pose problems for a jet-only model fit. 

These results further highlight the importance of obtaining simultaneous radio data 
with the smallest errors and/or as many simultaneous bands as possible to place tight 
contraints on the radio spectral slope. Sampling the optical/IR spectral region must now 
include the mid-IR and UV to sufficiently constrain theoretical models. Only then can we 
begin to have confidence in the location of the power-law break and how this varies over the 
duration of an outburst (i.e. with varying X-ray flux). 

Timing analysis during the soft (V-band) and hard states (//-band) show the variability 
of GX 339-4 is consistent with a simple power-law over timescales ~4-230 days. We see no 
evidence of characteristic timescales (such as power- law breaks or QPOs). Other studies 



that have obser yed QPOs during t he hard (iMotch et al 



1985 



Gandhi 



20101 ) and soft (IMotch et al. 



2009 : 



Gandhi et al 



19851 ) state have been for timescales much less than what 
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we've observed here (< 250 sees). In future work we will obtain data with higher temporal 
resolution in an attempt to bridge the timescale-gap between 250 sees and 4 days during 
both hard and soft X-ray states, and construct PSDs, cross-, and auto-correlation functions 
to fully explore variability timescales in GX 339-4 to see any evidence of characteristic 
timescale(s) related to accretion processes. 

We will continue to monitor GX 339-4 in the 0/IR with SMARTS for the foreseable 
future, including data at B-, V-, I-, J- H-, and ii'-bands, making it immediately public via 
the wetcl for use by the general astronomical community. 

We are grateful to our SMARTS ANDICAM observers, Juan Espinoza, Alberto 
Miranda, and David Gonzalez, and to Suzanne Tourtellotte for reducing the optical data. 
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Fig. 1. — V-band finding chart for GX 339-4. Field of view is 3.7' x 2.7'. 
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Fig. 3. — Photometric error of stars in the GX 339-4 field versus magnitude for optical 
and near- infrared bands. The solid line is a 4th-order polynomial fit. The y-axis is the 1-a 
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Fig. 4. — The SMARTS l.Sin optical and near-infrared light curve of GX 339-4 during 
2002-2010, inclusive. 
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Fig. 5. — 0/IR light curves of GX 339-4 during 2003 and 2005, highlighting the rebright- 
ening events. The solid line is a best-fit line to the data immediately before the rise of the 
rebrightening event. The slopes of these lines are given in Table HI 
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Fig. 6. — Hardness-intensity diagram of GX 339-4 for 2002-2010. Each separate outburst is 
shown in different colors. 
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Fig. 8. — i7-band light curve (top) and RXTE data. We categorized hard/soft state iJ-band 
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ratio is > 0.6 (hard state) and red solid triangles when the hardness ratio is < 0.6 (soft 
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Fig. 9. — Data are represented using the same criterion as Figure |H] and categorized into 
hard and soft states only for days when there were simultaneous RXTE data. 
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Hard-to-Soft Transitions (2002, 2010) 



Sott-lo-Hard Transitions (2003, 2005, 2007) 





Fig. 10. — Evolution of soft-to-liard and hard-to-soft state transitions compared to the V/ H 
correlation (open circles). Left : Hard-to-soft transition: Red line is 2002 data, green line 
is 2010 data. The soft-to-hard transition traces a unique path relative to the hard-to-soft 
transition. Right: Soft-to-hard transition: Red line is 2003 data, green line is 2005 data, 
and blue is 2007. In 2003 and 2005, GX 339-4 fades to the faintest end of the branch, but 
remains in the hard-state branch in 2007 (and 2008/2009). 
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Fig. 11. — Spectrophotometric standards (left) LTT 4816 and (right) EG 274. Triangles 
are, from left to right, SMARTS B, V, R, and / data calibrated using methods outhned 
in Sections 12.1.11 and 14.11 Vertical errorbars represent the 0.02 mag calibration error, and 
horizontal errorbars show the full-width-half-maximum of the filter bandpass. 
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Fig. 12.— Soft-state SEDs of GX 339-4 for V = 16 and 17 mag. /-, J-, and i7-band data 
were extrapolated from the power-law fits to the correlations shown in Figure [71 The errors 
are the combined photometric and reddening errors as outlined in Section 14.11 The dashed 
line represents the spectral slope expected for a viscous-heated accretion disk (F oc z/^/^). 
Dotted lines are the best-fit power-laws to the OIR data; ay=iQ = 1.8±0.1; at,=i7 = 1.6±0.1. 
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Fig. 13. — Hard state SEDs of GX 339-4 for V = 19.5 - 15.0 mag. Dashed line represents a 
viscous-heated accretion disk {F oc z/^/^). The sohd hne is a best-fit power-law to the 0/IR 
data {a = 1.3 ±0.1) aX V = 19.5 mag. This line is repeated as dotted lines and shifted up to 
match the V-band, highlighting that the V and / bands remain on the same spectral slope 
as GX 339-4 increases in flux. The J and H band spectral slope decreases with increasing 
brightness but the near-IR flux also increases suggesting an additional, non-thermal source 
becomes more dominant as GX 339-4 increases in flux. 
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Hard State, V - 15.73 mag 



Hard State, V - 15.14 mai 





log(Frequency) 
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Fig. 14. — Power-law model fits to the hard-state SED for V = 15.73 {left) and 15.14 
(right), when simultaneous 0/IR and radio data were obtained on UT 03 and 18 April 
2002, respectively. The dotted line is a power-law fit to the radio data only [left: a = 
0.11 ± 0.10; right a = 0.14 ± 0.01). The error on the radio data on 03 April 2002 is large 
enough such that the extrapolation may lie well above and below the OIR data. On 18 
April 2002, the radio spectrum lines up with the near-IR. The solid line is a power-law fit 
to the radio, J and if-bands [left: a = 0.18 ± 0.01; right a = 0.17 ± 0.01). In this case, 
the radio and near-IR line up on the same spectral slope (within errors) on both dates and, 
therefore, originate from the same fiux source, most probably jets. The V and I bands are 
more consistent with a thermal fiux source (as shown in Figure [T51) . 
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Fig. 15. — Power spectral densities of the hard state H-hand, and soft state ^-band, hght 
curves. The PSD of the observed data is shown by the sohd jagged hne, while the underlying 
power-law model is shown by the dashed straight line. Points with error bars correspond to 
the mean value of the PSD simulated from the underlying power-law model. The error bars 
are the standard deviation of the distribution of simulated PSDs. The broadband power 
spectral density is best described by a simple power law with a slope —1.8 at H-hand, and 
—2.0 at l^-band. See Table [6] for the PSD fit results for the \^-band during the hard state 
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Table 1. Optical comparison star magnitudes. Errors are 1-a. 

Number RA Dec V I {V - 1) 

(J2000.0) (J2000.0) (mag) (mag) (mag) 

1 17:02:52.0 -48:47:09.5 17.54 ± 0.06 15.26 ± 0.07 2.28 

2 17:02:49.9 -48:47:17.1 17.43 ± 0.05 15.49 ± 0.07 1.94 

3 17:02:51.8 -48:47:21.8 16.99 ± 0.05 15.62 ± 0.06 1.37 

Note. — Star numbers are in accordance with the optical finding chart (Figured]). 
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Table 2. Near-infrared comparison star magnitudes. Errors are 1-a. 

Number RA Dec J H 

(J2000.0) (J2000.0) (mag) (mag) 

1 17:02:52.0 -48:47:09.5 13.35 ± 0.09 12.49 ± 0.09 

2 17:02:49.9 -48:47:17.1 13.88 ± 0.09 13.02 ± 0.10 

Note. — Star numbers are in accordance with the IR finding chart (Figure [2]). Our 
magnitudes are consistent with the 2MASS catalog within our photometric errors. 



Table 3. Calibrated V-, I-, J-, and H-hand magnitudes, fluxes and their errors of GX 339-4 from 2002 - 2010, 
inclusive (abbreviated). The JD is that at the start of the observing sequence. The full version can be accessed online. 





year 


JD 


V mag 


V err 


Vo fiux 


Vo err 


/ mag 


/ err 


7o fiux 


/q err 


,7 mag 


,7 err 


Jo flux 


Jo err 


H mag 


H err 


Ho flux 


Ho err 




(-2450000.0) 


(mag) 


(mag) 


(Jy) 


(Jy) 


(mag) 


(mag) 


(■ly) 


(Jy) 


(mag) 


(mag) 


(Jy) 


(Jy) 


(mag) 


(mag) 


(Jy) 


(Jy) 


2002 


2297.85400 


. 


. 


. 


. 


17.730 


0.030 


0.00103 


0.00010 


. 
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. 


. 


16.080 


0.070 


0.00069 


0.00007 




2298.87070 
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16.640 
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0.00041 
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0.00006 


- 


- 






16.360 


0.090 
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0.00116 


0.00011 
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0.00007 
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0.00013 


- 


- 


- 
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Table 4. Slopes of best-fit lines to 2003 and 2005 light curves (see Figure E]) immediately 
before the rebrightening events. Errors stated are 1-cr. 

Year V I J H 

(mag/day) (mag/day) (mag/day) (mag/day) 

2003 0.010 ± 0.002 0.008 ± 0.002 0.007 ± 0.003 0.003 ± 0.002 
2005 0.012 ± 0.001 0.010 ± 0.001 0.009 ± 0.001 0.007 ± 0.001 
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Table 5. Power-law fit parameters to 0/IR correlations shown in Figure [71 Soft-state 

branch requires only a single power-law fit {Fy = C.F^^). The hard-state branch requires a 

broken power-law {Fy = C.F^^ for x < Xbreak, and Fy = C.F^\F"^~^^j^-^ for x > Xbreak)- 

Errors are l-o". 

Data/Branch ai Q!2 Xbreak C 

(mag) (mag) 

V/I hard 1.039 ± 0.001 0.890 ± 0.005 16.07 ± 0.03 -2.06 ± 0.02 

V/J hard 1.242 ± 0.002 0.858 ± 0.006 16.20 ± 0.01 -7.04 ± 0.04 

V/H hard 1.362 ± 0.002 0.853 ± 0.005 16.30 ± 0.01 -9.89 ± 0.04 

V/I soft 0.893 ± 0.002 n/a n/a 0.71 ± 0.03 

V/J soft 0.796 ± 0.004 n/a n/a 1.41 ± 0.06 

V/H soft 0.722 ± 0.003 n/a n/a 2.25 ± 0.06 
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Table 6. A summary of the best-fit power-law parameters to the power spectral densities 
of GX 339-4 during the hard and soft states. Errors are l-o". 

Filter /X-ray state Start JD End JD Slope (a) 

(-2450000.0) (-2450000.0) 



V soft state 


3224.70050 


3478.79434 


^•U_0.5 


V hard state 


4498.88162 


4746.54995 


1 8+0-3 


H soft state 


3224.70050 


3478.79434 


1 7+0.2 
-■-•'-0.6 


H hard state 


4498.88162 


4746.54995 


1 8+0-3 



